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1.  In  this  paper  we  conduct  a  novel  investigation  of 
tropical  cyclone  genesis  and  intensification  from  the 
perspective  of  helical  features  of  atmospheric  flows  of 
different  scales,  which  contribute  to  the  organization 
of  the  cyclone.  Using  the  data  from  a  near  cloud¬ 
resolving  numerical  simulation,  helical  characteristics 
of  the  velocity  field  are  calculated  and  analyzed  during 
a  process  of  formation  of  tropical  depression  vortex 
and  its  subsequent  intensification  up  to  a  mature  hur¬ 
ricane  stage.  An  evolution  of  large-scale  vortex  insta¬ 
bility  is  followed,  which  develops  on  the  background  of 
intense  cloud  moist  convection.  The  instability 
progresses  by  merging  of  small-scale  helical  convec¬ 
tive  structures  and  results  in  the  formation  of  larger 
and  more  intense  spiral  vortices.  Distributions  of 
helicity,  kinetic  energy  and  enstrophy  have  been 
obtained.  It  has  been  discovered  that  the  process  is 
accompanied  by  a  break  of  the  reflection  symmetry  of 
three  dimensional  moist  convective  turbulence  that  is 
manifested  in  the  generation  of  nonzero  mean  helicity. 
With  such  a  break  of  symmetry  a  theoretical  hypothe¬ 
sis  is  proposed  to  explain  the  self-organization  of 
atmospheric  processes  under  specific  features  of  heli¬ 
cal  turbulence,  which  provide  an  energy  transfer  from 
small-scales  structures  to  larger  vortices.  The  goal  of 
this  collaborative  Russian— American  investigation, 
which  started  recently,  is  a  clarification  of  the  role  of 
helicity  in  tropical  cyclogenesis  and  intensification. 

To  mathematically  describe  the  helical  features  of 
atmospheric  flows  it  is  useful  to  examine  a  well-known 
fluid-dynamical  characteristic  of  such  flows,  called 
helicity  of  the  velocity  field  [  1  ] .  Helicity  of  the  velocity 
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field  is  a  pseudoscalar  quantity  defined  as  the  dot 
product  of  velocity  and  vorticity  vectors  (see,  e.g., 
review  article  [2]  and  references  therein  on  helicity 
and  helical  turbulence).  A  volume  integral  calculated 
in  a  specific  space  domain 

H  =  jV-  rotVc/r  (1) 

gives  the  total  helicity  of  vortex  system.  A  non-vanish¬ 
ing  mean  helicity,  (H)  ^  0,  implies  the  symmetry  break 
of  turbulence  with  respect  to  coordinate  system  reflec¬ 
tions  [3].  Its  sign  determines  the  predominance  of  the 
left-handed  or  the  right-handed  spiral  motions.  For 
instance,  positive  mean  helicity  will  be  generated  in 
the  moist  atmosphere  under  the  predominance  of 
cyclonic  updrafts  and/or  anticyclonic  downdraft 
motions.  Similarly,  negative  helicity  will  be  generated 
for  the  case  of  anti-cyclonic  updrafts  and/or  cyclonic 
downdraft  flows. 

In  tropical  cyclone  investigations,  helical  features 
of  the  velocity  field  have  not  been  highlighted  until 
very  recently.  As  far  as  we  are  aware,  only  papers  [4—7] 
can  be  indicated  on  this  subject.  In  [4,  5]  calculations 
of  helicity  for  mature  hurricanes  were  carried  out  by 
regional  atmospheric  models  whose  space  resolution 
did  not  allow  resolving  the  most  energetic  scales  of 
cumulus  cloudiness,  while  in  [6,  7]  an  analysis  of 
helicity  values  was  performed  in  the  context  of  its  pos¬ 
sible  role  in  increasing  the  ability  of  developing  hurri¬ 
cane  to  resist  the  negative  impact  of  the  ambient  verti¬ 
cal  wind  shear.  In  papers  [6,  7]  helicity  was  calculated 
using  tropospheric-deep  dropsonde  soundings  carried 
out  by  reconnaissance  aircrafts  under  investigations  of 
eight  tropicai  cyclones  of  1998—2001  seasons  during 
wide-ranging  campaign  CAMEX  (Convection  and 
Moisture  Experiment)  [8]  organized  by  NASA 
(National  Aeronautics  and  Space  Administration). 

There  exists  a  great  amount  of  investigations  con¬ 
cerned  with  examining  the  role  of  cloud  convection 
(without  consideration  of  its  helical  features)  in  tropi¬ 
cal  cyclone  formation.  There  is  a  growing  consensus 
that  deep  cumulonimbus  convection  of  2—20  km  hor¬ 
izontal  scales,  which  transfers  sensible  and  latent  heat 
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from  the  underlying  surface  throughout  the  tropo¬ 
sphere  layer  in  the  tropics,  represents  the  main  mech¬ 
anism  to  intensify  a  pre-existing  cyclonic  circulation 
on  the  atmospheric  mesoscales  (-200  km)  to  a  vortex 
of  hurricane  strength  [9—11].  However,  at  present,  a 
clear  consensus  of  opinion  has  not  yet  emerged  con¬ 
cerning  a  scenario  of  such  transformation  and  physical 
mechanisms  contributing  to  it. 

Nearly  thirty  years  ago  in  paper  [12]  was  proposed 
a  mechanism  for  intensification  and  sustaining  of 
large-scale  vortex  disturbances  in  the  atmosphere  due 
to  energy  transfer  from  small-scale  helical  convective 
turbulence — the  so  called  turbulent  vortex  dynamo. 
Only  recently  has  it  become  possible  to  test  this 
hypothesis  with  a  physically  consistent  data  set.  And 
only  recently  have  studies  been  conducted  to  examine 
reasonably  high  horizontal  resolution  (-1—3  km  of 
horizontal  scale  and  less)  numerical  simulations  of 
tropical  cyclone  formation  that  possess  an  adequate 
representation  of  both  the  deep  cumulus  and  strati¬ 
form  stages.  Between  these  works  [13]  offered  a  new 
scenario  of  hurricane  formation  based  on  self-organi¬ 
zation  of  convective  processes  in  an  otherwise  favor¬ 
able  tropical  environment.  Although  helical  features 
of  these  simulated  flows  were  not  taken  into  consider¬ 
ation  within  the  framework  of  paper  [13],  self-organi¬ 
zation  of  vortical  (!)  convection  was  observed  similar 
to  “helical”  scenario  [2],  namely,  as  an  enlargement  of 
vortex  structures  from  the  size  of  individual  rotating 
cumulus  clouds  in  the  model,  their  induced  concen¬ 
tration  of  absolute  angular  momentum  on  the  system 
scale  circulation,  and  their  merging  with  each  other  to 
yield  newly  forming  larger  vortices  and  an  intensifying 
circulation  on  the  system  scale. 

In  this  work  we  use  the  velocity  fields  obtained  in 
[13]  to  calculate  and  analyze  helical  characteristics  of 
the  cyclogenesis  and  intensification  process  for  the 
problem  as  posed  by  [  1 3] . 

2.  A  numerical  approach  for  investigation  of  large- 
scale  helical  instability  in  the  atmosphere  was  pro¬ 
posed  in  [14,  15]  and  first  applied  in  a  simpler  case, 
namely,  to  simulate  helical-vortex  effects  in  laminar 
Rayleigh— Benard  convection  by  use  of  an  additional 
helical  force.  Simulation  [14,  15]  demonstrated  quali¬ 
tatively  new  effects  which  appeared  in  a  crisis  manner. 
A  nonzero  mean  helicity  ( H)  ^  0  of  the  flow  was  gen¬ 
erated,  and  after  exceeding  some  its  value,  the  large- 
scale  vortex  instability  appeared.  A  key  role  in  the  pro¬ 
cess  the  vertical  component  of  velocity  played,  and  the 
corresponding  evolutionary  equation  in  mathematical 
model  [14,  15]  closed  a  positive  feedback  loop.  The 
instability  evolution  observed  as  an  enlargement  of 
horizontal  scales  of  structures  what  happened  by 
merging  of  helical  vortex  cells.  The  process  was 
accompanied  by  flow  intensification  in  newly  forming 
larger  vortices,  similarly  to  that  observed  in  atmo¬ 
spheric  modeling  [13],  which  was  based  on  hydro¬ 
thermodynamic  equations  for  the  atmosphere. 


It  should  be  particularly  pointed  out  about  a  dis¬ 
tinct  enhancement  of  the  heat  transfer  discovered  in 
experiments  [14,  15].  The  heat  flux  through  a  layer 
increased  with  an  increase  in  mean  helicity  of  the  flow. 
A  sharp  intensification  of  heat  transfer  observed  after 
exceeding  the  threshold  of  stability,  and  later  during 
the  instability  evolution,  after  each  merging  of  vortex 
structures.  Thus,  the  effectiveness  of  heat  transfer 
increased  due  to  two  factors,  namely,  increase  in  the 
mean  helicity  and  decrease  in  a  cell  number  of  the 
convective  system.  Such  energetic  expedience  sug¬ 
gests  to  us  that  a  most  effective  “removal”  of  accu¬ 
mulating  heat  might  be  possibly  one  of  the  reasons 
for  “helical”  self-organization  of  convection  in  trop¬ 
ical  cyclogenesis. 

3.  In  the  present  paper,  ideas  and  methods  pro¬ 
posed  previously  for  calculations  and  analysis  of  heli¬ 
cal  characteristics  [14,  15]  are  applied  to  post-process¬ 
ing  of  velocity  fields  obtained  in  [  1 3]  by  use  of  numer¬ 
ical  meteorological  model  RAMS  (Regional 
Atmospheric  Modeling  System)  developed  at  Colo¬ 
rado  State  University.  Detailed  information  on  the 
governing  equations  of  hydro-thermodynamics  of  the 
atmosphere,  parameterizations  of  turbulent  processes, 
specific  characteristics  of  the  RAMS  model  configura¬ 
tion,  initial  and  boundary  conditions  is  given  in  paper 
[13]  as  well  as  references  to  corresponding  works. 

3. 1.  Let  us  note  further  a  feature  that  is  of  crucial 
importance  for  the  problem  under  consideration.  The 
velocity  fields  used  for  calculations  of  helical  and  inte¬ 
gral  characteristics  in  this  work  were  obtained  in  [13] 
by  use  of  three-dimensional  nonhydrostatic  numerical 
modeling  system  comprising  time -dependent  equa¬ 
tions  for  all  three  components  of  velocity  (including 
the  vertical  one,  see  item  2  on  this  subject),  pressure, 
potential  temperature,  total  water  mixing  ratio,  and 
cloud  microphysics.  RAMS  utilizes  an  interactive 
multiple  nested  grid  scheme  which  allows  explicit  rep¬ 
resentation  of  cloud-scale  features  within  the  finest 
grid  while  enabling  a  large  domain  size  to  be  used, 
thereby  minimizing  the  impact  of  lateral  boundary 
conditions.  For  all  numerical  experiments  [13]  three 
nested  grids  were  used.  A  standard  radiation  boundary 
condition  was  used  at  the  lateral  boundaries,  which 
assumes  that  disturbances  reaching  the  boundaries 
move  as  linearly  propagating  gravity  waves.  A  standard 
Rayleigh  friction  layer  was  included  at  upper  levels  in 
order  to  minimize  reflection  of  gravity  waves  from  the 
top  of  the  model.  All  microphysical,  radiative,  and  dif¬ 
fusive  parameters  were  the  standard  ones  employed  for 
tropical  summer  conditions  [13].  The  initial  tempera¬ 
ture  distribution  in  [13]  was  the  mean  Atlantic  hurri¬ 
cane  season  sounding  which  was  representative  of  so- 
called  “non-Saharan-air-layer”  air. 

3.2.  We  performed  calculations  and  analysis  of  heli¬ 
cal  and  integral  characteristics  for  six  of  nineteen 
numerical  experiments  [13],  which  are  presented  in 
the  table.  For  the  purpose  of  expedience,  postprocess¬ 
ing  of  the  model  data  was  carried  out  on  the  finest 
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RAMS  numerical  experiments  [13J  analyzed  in  this  work 


No. 

Name 

of  experiment 

Max  v(m  s  *) 
at  z  =  4  km 

Description  of  experiment 

Al 

Control 

6.6 

Ax  =  Ay  =  2  km,  SST1  =  29°C.  Metamorphosis  to  surface  vortex  successful. 
Becomes  miniature  tropical  cyclone  by  approximately  60  h.  Mean  near-surface 
tangential  wind  ~12  m  s_1  at  24  h,  and  46  m  s_l  at  72  h. 

A2 

3  km 

6.6 

Ax  =  Ay  =  3  km,  SST  =  29°C.  Metamorphosis  to  surface  vortex  successful. 

Mean  near-surface  tangential  wind  ~13  m  s_1  at  24  h,  and  46  m  s_1  at  72  h. 

B3 

CAPE-less2 
(3  km) 

6.6 

Ax  =  Ay  =  3  km,  SST  =  29°C.  Low-level  moisture  decreased  by  2  g  kg-1. 
Metamorphosis  successful,  but  slower  rate  of  development.  Mean  near-surface 
tangential  wind  ~9  m  s_1  at  48  h. 

Cl 

No  vortex 

- 

Ax  =  Ay  =  3km,  SST  =  29°C.  No  initial  vortex.  Nosurface  development  whatsoever. 

C3 

Weak  vortex 

5.0 

Ax  =  Ay  =  3  km,  SST  =  29°C.  Metamorphosis  successful,  but  slower  rate  of  devel¬ 
opment.  Mean  near-surface  tangential  wind  ~9  m  s_1  at  72  h.  Circulation  very 
asymmetric  even  at  72  h. 

El 

Zero  Coriolis 

6.6 

Ax  =  Ay  =  3  km,  SST  =  29°C.  Coriolis  parameter  set  to  zero  (f=  0).  Metamorpho¬ 
sis  successful.  Develop  surface-concentrated  vortex  as  in  Al,  but  no  subsequent 
intensification  observed  through  72  h. 

Note:  1  SST — Sea  Surface  Temperature. 

2  Convective  Available  Potential  Energy. 


computational  grid  for  subsequent  times  with  a  time 
increment  of  10  min  during  72  h  of  numerical  experi¬ 
ment.  Characteristics  were  calculated  in  the  computa¬ 
tional  domain  of  276  x  276  x  20  km  in  Cartesian  coor¬ 
dinates  (x,  y,  z)  by  use  of  uniform  finite-difference  grid 
with  increments  Ax,  Ay,  A z-  In  all  six  experiments  the 
vertical  increment  in  the  post  analysis  was  identical 
and  equal  to  500  m;  the  horizontal  increments  were 
Ax  =  Ay  =  3  km,  except  experiment  A1  with  Ax  =  Ay  = 
2  km. 

3.3.  To  analyze  the  process  of  self-organization  of 
moist  atmospheric  convection  observed  under  condi¬ 
tions  of  tropical  cyclogenesis  as  posed  by  [  1 3] ,  a  set  of 
helical  characteristics  was  computed,  as  well  as  some 
other  integral  characteristics  of  the  velocity  field 
which  were  applied  in  [  14,  15]. 

The  following  characteristics  were  calculated  with 
a  time  increment  10  min  during  the  whole  72-hours 
evolution  of  tropical  cyclone  to  obtain  results  dis¬ 
cussed  in  this  paper:  three-dimensional  relative  helic- 
ity  density 

Huk  =  (  V  ■  rotV),  /  4,  (2) 

two  other  important  characteristics  of  the  velocity 
held — enstrophy  and  kinetic  energy  densities: 


1  2 

Hj.k  =  -(rofV),^*, 


Euk  =  \<y)u,k,  (3) 


as  well  as  mean  (volume-integrated)  values  of  helicity, 
enstrophy  and  kinetic  energy  integrated  over  the  whole 
computational  domain  276  x  276  x  20  km  and  normal¬ 
ized  by  number  of  grid  points: 

(H),  <8>,  (E).  (4) 

4.  Before  we  begin  the  discussion  of  results,  it  is 
helpful  to  remember  the  characteristic  values  of  wind 


velocity  for  three  stages  of  tropical  cyclone  strength  so 
we  may  readily  compare  them  with  data  in  the  table.  At 
the  stage  of  tropical  depression  the  velocity  of  near 
surface  tangential  wind  does  not  exceed  17  m  s-1, 
velocity  values  within  17— 33  m  s_l  correspond  to  the 
tropical  storm,  more  than  33  m  s-1  to  the  hurricane 
strength. 

In  numerical  experiments  Al,  A2,  B3,  C3,  and  El 
from  [13]  a  transformation  of  an  initial  mid-tropo- 
spheric  vortex  disturbance  into  a  surface-concentrated 
tropical  depression  vortex  was  observed.  Some  unfore¬ 
seen  differences  in  both  duration  of  the  process  and 
intensity  of  the  formed  vortex  were  found  between 
these  experiments.  However,  experiment  El,  in  the 
absence  of  the  Coriolis  force,  did  not  reveal  any  fur¬ 
ther  intensification  of  the  surface-concentrated  vortex 
during  72  h.  In  case  Cl,  without  an  initial  vortex  dis¬ 
turbance,  no  system-scale  surface  development  was 
observed. 

In  experiments  Al ,  A2,  B3,  and  C3  an  influence  of 
different  factors  (horizontal  resolution,  convective 
process,  intensity  of  initial  vortex)  on  cyclogenesis  was 
studied,  and  similar  tendencies  were  observed  in  evo¬ 
lution  of  helicity,  enstrophy  and  kinetic  energy. 

Let  us  now  discuss  the  results  of  numerical  experi¬ 
ment  A2  “3  km.” 

Under  examination  of  evolution  of  three  dimen¬ 
sional  helicity  field  (2),  vertical  velocity  and  vertical 
vorticity  were  analyzed  also .  This  study  readily  allowed 
an  identification  of  the  formation  of  convective  struc¬ 
tures  and  determining  their  rotational  signature,  i.e., 
cyclonic  or  anticyclonic. 

The  undertaken  modeling  study  showed  the  exist¬ 
ence  of  intense  helical  convection,  whose  ascending 
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Fig.  1.  Evolution  of  kinetic  energy  and  helicity  integrated  over  the  computational  domain  and  normalized  by  number  of  grid 
points:  experiment  A2  (left  column — (a,  b)),  experiment  El  (right  column — (c,  d)). 


vertical  velocities  sometimes  exceeded  30  m  s-1  during 
a  portion  of  the  numerical  experiment,  starting  near 
t=  -10—15  h  [13].  We  observed  the  predominance  of 
convective  vortex  structures  with  cyclonic  rotation  of 
ascending  flow  (i.e.,  positive  helicity).  A  spectmm  of 
convective  motions  was  clearly  observed,  whose  horizon¬ 
tal  scales  changed  from  4—9  up  to  15—30  km,  while  local 
helicity  values  in  them  were  0.002—0.004  m  s-2.  The 
most  intense  of  these  structures  reached  through  the 
bulk  of  the  troposphere.  Hot  vortical  plumes  that  rise 
through  the  bulk  of  the  troposphere  have  been  called 
“vortical  hot  towers.”  Shortly  after,  a  process  of  con¬ 
vective  vortices  merging  started,  and  structures  whose 
horizontal  sizes  reached  10—30  km  were  dominated, 
their  helicity  exceeded  0.008  m  s-2.  During  the  process 
of  merging  one  could  observe  an  increase  in  the  back¬ 
ground  helicity  in  adjacent  areas. 

As  a  result  of  further  subsequent  mergers,  a  surface - 
concentrated  vortex  of  tropical  depression  strength 
(with  mean  tangential  wind  about  13  m  s_l)  formed 
after  24  h  of  simulation.  During  further  development 
into  a  hurricane  power  vortex  (46  m  s-1  at  t  =  72  h) 


even  more  intense  vortical  hot  towers  appeared,  they  pos¬ 
sessed  helicity  values  of  0.5— 1.0  m  s-2. 

Evolution  of  kinetic  energy  and  helicity,  calculated 
by  use  of  (4),  during  the  process  of  hurricane  forma¬ 
tion  is  shown  in  Figs,  la,  lb.  Enstrophy  behavior  is 
qualitatively  similar  to  that  observed  for  kinetic  energy 
in  Fig.  la. 

All  pictures  distinctly  show  an  initial  period  of 
approximately  10  h,  which  is  needed  for  development 
of  intense  small-scale  helical  convection  and  starting 
of  the  process  of  merging  of  convective  structures. 
During  next  1 0—  1 5  h,  a  self-organization  process  con¬ 
tinues  and  leads  to  formation  of  tropical  depression 
vortex  which  has  essentially  larger  scales.  A  local  peak 
in  all  pictures  corresponds  to  this  event  at  approxi¬ 
mately  24—25  h.  Further  intensification  of  tropical 
depression  up  to  hurricane  power  till  the  end  of  simu¬ 
lation  time  at  72  h  can  be  traced  as  a  strong  increase  in 
kinetic  energy  and  helicity  values. 

It  should  be  pointed  out  that  between  the  three 
integral  characteristics  of  the  velocity  field,  only  helic¬ 
ity  demonstrates  qualitatively  different  behavior.  Dur¬ 
ing  the  first  15—17  h  the  integral  helicity  is  close  to 
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zero,  and  there  exist  short  time  periods  when  it 
becomes  negative.  After  approximately  1 8  h  of  simula¬ 
tion  time  helicity  becomes  distinctly  positive  and 
increasing.  Analysis  of  vertical  velocity  and  vertical 
vorticity  for  these  conditions  shows  that  such  positive 
helicity  is  a  result  of  the  predominance  of  locally 
cyclonic  updrafts.  Non-zero  helicity  means  a  break  of 
the  mirror  symmetry  of  atmospheric  turbulence  what 
may  result  in  generation  of  a  large-scale  instability  [3] . 

It  is  instructive  to  compare  the  kinetic  energy  and 
helicity  evolution  in  this  experiment  A2  (Figs,  la,  lb) 
against  that  of  experiment  E 1 ,  the  absence  of  the  Cori¬ 
olis  force  (Figs.  lc,  Id).  In  the  latter  case  when  the  sur¬ 
face-concentrated  vortex  forms,  there  is  no  further 
intensification.  As  we  can  see,  in  scenario  El  the 
helicity  behavior  is  cardinally  different:  the  oscilla¬ 
tions  are  of  significant  amplitude,  and  are  even 
accompanied  by  a  change  of  sign  when  the  largest  pos¬ 
itive  and  negative  values  are  comparable  to  their  abso¬ 
lute  values. 

It  is  important  to  point  out  that  the  values  of  posi¬ 
tive  helicity  generated  in  experiment  A2,  which 
resulted  in  hurricane  formation,  are  more  than  two 
times  higher  than  positive  helicity  values  in  case  E 1 . 
Furthermore,  no  change  in  the  sign  in  helicity  is 
observed  in  A2  after  the  tropical  depression  stage  is 
attained. 

The  results  presented  portray  an  important  signa¬ 
ture  of  integrated  helicity  in  the  evolution  of  large- 
scale  vortex  instability.  They  clearly  motivate  further 
investigations  of  spectral  characteristics  of  helical 
moist  convective  atmospheric  turbulence  in  condi¬ 
tions  of  tropical  cyclogenesis. 

The  practical  significance  of  these  investigations  is 
connected  with  a  possibility  to  use  helicity  as  an  indi¬ 
cator  of  large-scale  vortex  instability. 
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